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A REAL-GAS STUDY OF LOW-DENSITY WEDGE-INDUCED 


LAMINAR SEPARATION ON A HIGHLY COOLED 
BLUNT FLAT PLATE AT M*, = 12 

By John B. Anders and C. L. W. Edwards 
Langley Research Center 

SUMMARY 

A study has been made of real-gas wedge -induced laminar separation. A series 
of experimental tests on a blunt flat plate with a trailing-edge flap was made at a free- 
stream Mach number of 12 and a free-stream Reynolds number of 10 4 /ft (3.3 X 10 4 /m). 
The tests were conducted for stagnation enthalpies ranging from 1465 Btu/lbm 
(3.41 MJ/kg) to 2030 Btu/lbm (4.73 MJ/kg). A calculation technique is developed to 
predict the chordwise extent of wedge -induced laminar separation and agrees reasonably 
well with the present real-gas experimental data. The computed effect of free-stream 
Reynolds number and total enthalpy on the extent of separation is demonstrated. Sample 
calculations are shown for equilibrium and frozen flow for a range of total enthalpies. 

INTRODUCTION 

The occurrence of laminar -boundary-layer separation on maneuverable reentry 
vehicles results in large changes in aerodynamic characteristics, heat transfer, and con- 
trol forces. This problem has been studied by many investigators and much theoretical 
and experimental information exists on laminar separated flows. (See refs. 1 to 22.) 
However, all the theoretical investigations, and all the experimental investigations with 
the exception of Kuehn's work (refs. 23 and 24) have considered the fluid to be a perfect 
gas. Also, with the exception of Rogers and Berry (ref. 25) and Kuehn (refs. 23 and 24), 
only relatively high-density flows have been investigated, usually with local Reynolds 
numbers above 10^. The higher kinetic energies of reentering vehicles are normally 
associated with the early low -density region of the flight path. Therefore, the regions 
of low density and high enthalpy are linked together for this application. 

The extension of present theories to low-density real-gas conditions invariably 
seems to violate assumptions made in their development. Similarly, the extrapolation of 
experimental results to low-density real-gas conditions is, perhaps, unfounded without 
some experimental investigation in this region. The present study was begun to check 



experimentally certain semiempirical perfect -gas relations in the real -gas region and 
to see whether certain gross predictions of low -density real-gas separation could be 
made by utilizing simple physical concepts. Such a study would be useful in demon- 
strating trends as well as pointing the way to future research. 

The model chosen for this study was a blunt-nosed wing -flap combination. All 
experimental testing was done in the Langley 1-foot hypersonic arc tunnel at a nominal 
Mach number of 12 and a nominal free-stream Reynolds number per foot of 1.4 x 10^ 
(4.59 x 10^ per meter). The Knudsen number based on nose diameter was approximately 
0.06. 


SYMBOLS 


A,B 

constants in equations (6) and (7) 

a 

speed of sound 

n 

ii 

3 HI 



constants in equations (11), (12), and (13) 

C P>P 

plateau pressure coefficient, ^ 

y o 

c f,o 

local-skin-friction coefficient 

d 

model nose diameter 

G - 1.7208 

^(^,0.3859) 

Ho 

form factor 

h 

enthalpy 

h t,2 

stagnation enthalpy at model nose 

k 

nose drag coefficient 

k l 

constant in equation (19) 

L 

model length to wedge hinge line 
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length of dividing streamline 


length over which pressure rise to plateau takes place 


Mach number 


defined by equations (14a) and (14b) 


static pressure 


bluntness induced pressure 


stagnation pressure on model nose 


viscous induced pressure increment 


dynamic pressure 


specific gas constant 


Reynolds number per foot (per 0.3048 meter) based on free -stream 
p V 

conditions, °° °° 

M oo 


free -stream Reynolds number based on plate length, 


PooVooL 


free -stream Reynolds number based on longitudinal distance 


Reynolds number based on conditions at beginning of interaction, 


distance variable from reference 9 


temperature 

velocity ratio along dividing streamline, — 

v 2 

velocity component 


longitudinal coordinate (see fig. 2) 


3 



z real-gas compressibility factor 

0 = 

y ratio of specific heat 

6 angle of dividing streamline with plate surface (see fig. 10) 

6* displacement thickness 

e function of y, e ~ 1 — 0.Q048 

(y - l ) 2 

6 wedge deflection angle 

ix coefficient of viscosity 

£ defined by equation (17) 

correction factor in equations (4) and (5) 
p mass density of air 

c p defined by equation (18) 

X viscous interaction parameter, 

Subscripts: 

aw adiabatic wall 

c computed value 

o condition at beginning of interaction 

p condition in plateau region 

r condition at reattachment 

s condition at separation 
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t 


stagnation condition 


w condition on model wall 

°o free-stream condition 

2 condition at edge of boundary layer 

APPARATUS 

The electric -arc -heated gas in the tunnel is expanded in a 5°-half-angle, conical 
nozzle to a nominal Mach number of 12. The inviscid core of this facility has been sur- 
veyed and was found to be approximately 5 inches (12.7 cm) in diameter. Test durations 
as long as 15 minutes have been utilized in this facility. A detailed description of the 
facility is presented in reference 26. 


Model 

The model used in the tests was a water-cooled, hemicylindrically blunted flat 
plate with a variable -angle trailing-edge flap. The plate had a nose radius of 1/8 inch 
(0.32 cm) and a chord of 2^ inches (6.35 cm) and completely spanned the test section. 

The flap had a chord of 1^ inches (3.81 cm) and a span of 4 inches (10.16 cm). The 
model was constructed of brass with a copper leading edge. The flap was designed for 
an angular variation of 20° to 45° and also for 0°. The model was also fitted with 
uncooled, detachable, stainless -steel end plates mounted 4 inches (10.16 cm) apart. 

A photograph of the model mounted in the test section is presented in figure 1. Figure 2 
shows the end plates on the model and the position of the instrumentation. The model 
was instrumented with 4 chromel-alumel thermocouples and 31 pressure orifices. The 
pressure orifices were connected to ionization-type pressure sensors and both the 
thermocouple and ionization gage outputs were recorded on oscillograph film recorders. 

Test Conditions 

The present series of tests were conducted in air for a nominal Mach number 
of 12, a stagnation enthalpy range from 1465 Btu/lbm (3.41 MJ/kg) to 2030 Btu/lbm 
(4.73 MJ/kg), a free-stream Reynolds number per foot from 1.1 x 10^ to 1.9 x 10^ 

(3.6 x 10^ per meter to 6.2 x 10^ per meter), and a flap angle range from 27° to 32°. 
Tunnel blockage prevented flap angles larger than 32°. Precise values for tunnel con- 
ditions during specific tests are presented in table I. Ratios of wall temperature to 
total temperature varied from 0.11 for the lowest enthalpy test to 0.09 for the highest 
enthalpy test. 
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Figure L- Model installed in 1-foot test section. End plates removed. L-67-8702 


TABLE I. - TUNNEL TEST CONDITIONS 


Btu/lbm MJ/kg atm 


Pt 

Poo 

Moo 

N/m 2 

atm 

N/m 2 

2.26 x 106 

5.47 x 10- 5 

5.54 

12.4 

2.06 

5.73 

5.80 

11.8 

2.25 

5.09 

5.16 

12.7 

2.21 

4.21 

4.26 

13.1 

2.21 

5.58 

5.65 

12.2 

2.17 

6.34 

6.42 

12.1 

2.14 

5.96 

6.04 

11.9 

2.13 

5.19 

5.26 

12.5 

2.16 

5.62 

5.70 

12.0 

2.07 

4.50 

4.56 

12.7 

1.89 

5.84 

5.92 

11.9 

2.40 

5.67 

5.74 

12.4 







Accuracy 


The maximum error in determining the total enthalpy has been found to be about 
±8 percent as reported in references 27 and 28. An orifice effect on pressure measure- 
ment has been investigated and reported in reference 29. This effect can cause a maxi- 
mum error in pressure measurements of -6 percent for the lowest pressures measured 
in these tests. Maximum error due to the orifice effect on higher pressures such as 
those at separation, plateau, or reattachment is considerably less than -6 percent (usually 
less than 1 percent). A maximum inaccuracy in instrumentation and reading of data has 
been found to be ±5 percent. The total uncertainty for the lowest measured pressure is 
then 5 percent and -11 percent. This uncertainty in pressure causes an uncertainty in 
the calculation of local properties along the plate. Cases were calculated by using the 
maximum error in pressure and this procedure was found to result in a 2-percent devia- 
tion in Cf 0 and a 4-percent deviation in local Mach number. These deviations were 
considered to be insignificant in the overall results. 

RESULTS AND DISCUSSION 

Blunt-Plate Pressure Distribution Without Separation 

Several tests were made with a 0° flap angle in order to determine the effect of the 
blunt nose on the pressure distribution over the plate. Figure 3 shows a typical measured 
pressure distribution with a 0° flap angle. No pressure measurements were made on the 
flap itself. Two theoretical distributions (refs. 30 and 31) are also shown. The method 
of Bertram and Blackstock (ref. 30) was found to give the best overall agreement with 
certain restrictions on y. The method itself is not specifically valid for a variable y 
flow; therefore, an attempt was made to represent the flow by some characteristic y. 
Kuehn (ref. 32) discusses the use of a characteristic y for blunt-nosed cylinders. He 
states that Wick of the NASA Ames Research Center has suggested that the value of y 
at the stagnation region should be used since the blast analogy describes the flow in terms 
of a blast wave originating at the nose. In addition to the inviscid blast wave pressure, 
the viscous induced pressure must be taken into account; and intuitively, this viscous 
induced pressure might be characterized by a locally constant y. From the analysis of 
reference 30, the inviscid and viscous pressures are given by the following expressions 
for Ap y small compared with p^ 

r -i2/3 

0.187e fy(y - 1)M?2^ + 0.74 (1) 

Poo x/d 
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and 


^Py _ 5 y^^-oo (2 

’ p ” " s /v p ~ 

where y in equation (1) is now evaluated at the model stagnation point and y in equa- 
tion (2) is a local value evaluated at pj^ and the stagnation -point entropy. The plate 
surface pressure is then given by 


p w + Ap v 


(3) 


These equations are used, without further qualification, for calculations presented later 
in this paper. 


30 


25 


20 


P/P=o 15 


10 


5 


Ref. 30 (eqs. (1) and (2) in text with y = 1.4) 

Ref. 30 (eqs. (1) and (2) in text with y in 

eq. (1) evaluated at model nose stagnation 
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Figure 3.- Blunt-plate pressure distribution. 

Separated Flow Results 

Test data were obtained by two methods: oil-flow studies and pressure distribu- 
tions. The partial-span flap of the present model leads to large three-dimensional effects 
in the separation bubble according to reference 33, and end plates were added to the model 
to minimize this spanwise bleed. A typical result from the addition of end plates is shown 
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in figure 4. The end plates significantly increased the chordwise extent of separation. 
With the exception of figure 4, all the data presented were taken with the end plates on 
and the effect of the end plates themselves is discussed in the light of the oil-flow pat- 
terns and spanwise pressure measurements. 



Figure 4.- End-plate effect. 


Oil-flow patterns .- The photographs of the oil-flow patterns were not of sufficient 
quality to warrant reproduction. However, figure 5 shows the pertinent features indicated 
by the original prints. Prior to a test, the surface of the plate was cleaned and dotted 
with drops of thin oil tinted with Prussian blue. The droplets were then observed during 
tunnel starting, approximately 2 minutes of testing time, and tunnel shutdown. The oil 
drops were observed to travel slowly back to a point on the plate surface, then to collapse 
and run together, and thus wet the plate surface. The complete process occurred in 3 to 
4 seconds and no change in the pattern could be detected during the remaining test time 
or during tunnel shutdown. The forward edge of the wetted area was assumed to be the 
separation line and the position of this line is approximately halfway up the pressure rise 
to plateau pressure. However, on two of the oil-flow tests, the separated region was very 
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small and the pressure distributions 
did not exhibit a definite plateau. 

For these tests the oil-flow pattern 
indicated that separation occurred 
nearly three -fourths of the way up 
the pressure rise to the inflection 
point. 

Spanwise pressure gradients. - 
All the oil-flow tests were made 
with end plates on, and from fig- 
ure 5 the flow appears to be 
approximately two dimensional 
near the model center line. The 
oil streaks on either end of the 
model near the end plates curve 
inward slightly and indicate some 
end-plate effects in this region. 
Outboard pressure measure- 
ments were made to check for 



Figure 5.- Sketch of typical oil-flow pattern. 


any spanwise pressure gradients in the separated region not indicated by the oil streaks. 
Unfortunately, one outboard orifice developed a water leak at the onset of the test program 
and had to be plugged. The remaining outboard orifice indicated center -line pressure 
within the reading accuracy in 6 of the 11 tests. The other 5 tests indicated an outboard 
pressure 22 percent to 27 percent lower than the center line. This result might be 
expected since the separation pressure rise for most of the tests occurred near the out- 
board orifice. A slight spanwise variation in separation length could cause a fairly large 
variation in the pressure level sensed by the outboard orifice. 


Cho rdwise pressure distributions. - Several typical pressure distributions are shown 
in figure 6, and table II lists the measured pressure distributions for all the tests. Fig- 
ure 6 shows the effect of flap angle and also the nose bluntness effects on the pressures 
near the leading edge. Also shown is the calculated inviscid wedge pressure for each 
test. The short chord of the flap raises the question of whether reattachment is being 
affected by flap chord. Since the measured maximum flap pressure is very near the cal- 
culated inviscid values, it is assumed that the flap chord length is sufficient to insure a 
reattachment process free from flap chord effects. Also, calculations based on the tech- 
nique presented later in this report indicate that reattachment should have occurred near 
the middle of the deflected flap for all except the 32° flap angle test. For the 32° test, 
reattachment appeared to occur at the top of the flap. In this instance the reattachment 
process was probably altered somewhat by the short flap chord. 
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Figure 6.- Ratio of pressure to tunnel free-stream static pressure. 


The pressure distributions were further reduced by expressing as a ratio the mea- 
sured pressure to the pressure on the plate with a 0° flap. Where 0° flap pressure dis- 
tributions were not available, extrapolated distributions obtained from the separated 
results were used. This type of data reduction allows a more meaningful comparison 
between tests as well as more clearly indicating the beginning of interaction. A typical 
plot of this type is shown in figure 7. The data for the 27° and 28° flap angles appear to 
show separation moving downstream with increasing flap angle. However, the 27° flap 
angle test had a higher Reynolds number and Mach number and a lower total enthalpy 
level and perhaps this fact accounts for the difference in the two tests. Nearly constant 
plateau regions were obtained in all the tests except those exhibiting only small separated 
regions. Rogers and Berry (ref. 25), in contrast to this result, did not find any evidence 
of a constant-pressure region in their investigation of low Reynolds number flow over a 
flat plate with forward -facing steps. The local Reynolds number on the plate for their 
tests ranged from 100 to 900 and indicated only a slight kink in the pressure distribution 
where plateau should have occurred. Gray (ref. 34) also questions the existence of a 
plateau in purely laminar flow. He indicates that this plateau is possibly a phenomenon 
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TABLE H.- MEASURED PRESSURE DISTRIBUTIONS 
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i 1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

0 

209.0 

193.1 

218.5 

234.4 

204.0 

199.0 

195.0 

213.0 

198.0 

219.1 

192.6 

209.3 

1.66 

13.58 

10.80 

14.74 

17.50 

20.17 

11.62 

11.59 

12.95 

12.18 

14.61 

14.98 

13.30 

2.41 

21.75 

14.37 

15.50 

17.80 

13.33 

11.84 

11.15 

12.31 

11.36 

13.75 

13.64 

12.20 

3.12 

9.85 

9.20 

10.33 

12.19 

9.55 

8.61 

8.50 

9.90 

8.90 

10.09 

9.57 

10.00 

3.75 

8.77 

8.05 

9.04 

10.62 

8.26 

7.26 

7.50 

10.91 

7.96 

9.06 

9.12 

8.40 

4.25 

8.17 

7.13 

8.27 

10.00 

7.55 

6.74 

6.74 

7.62 

7.02 

8.19 

8.00 

7.50 

4.75 

7.94 

7.36 

7.88 

10.62 

7.43 

6.33 

6.84 

9.14 

7.15 

5.85 

8.11 

7.60 

5.25 













5.75 

7.81 

7.13 

8.78 

16.25 

5.56 

6.43 

6.62 

7.74 

7.26 


8.44 


6.25 

8.78 

7.13 

10.33 

20.62 

7.43 

6.02 

6.62 

7.62 

6.67 

7.60 

7.55 

6.80 

6.75 

7.94 

7.36 

16.02 

25.00 

8.60 


7.50 

8.76 

7.15 

8.77 

10.35 

7.10 

7.25 

7.69 

7.36 

17.05 

27.50 

9.55 

13.50 

8.72 

10.28 

7.02 

8.77 

10.02 

6.70 

7.75 

11.06 

8.28 


29.70 

12.85 

7.26 

12.03 

12.81 

8.32 

11.40 

10.93 

6.60 

8.25 

14.30 

11.50 

22.23 

30.61 

16.15 

11.00 

15.12 

15.86 

10.89 

14.48 

12.74 

6.60 

8.75 

17.43 

14.70 

24.30 

37.20 

18.51 

11.63 

16.90 


14.51 

18.71 

17.11 

6.50 

9.25 

18.63 

15.29 

23.90 

34.40 

19.70 

16.39 

18.21 

18.91 

15.58 

19.88 

17.22 

6.00 

9.75 

20.20 

17.02 

24.55 

31.89 

20.75 

15.77 

19.42 

20.05 

17.10 

22.21 

18.36 

6.40 

10.50 

22.49 

20.10 

26.00 

33.75 

22.75 

18.06 

21.30 

22.10 

19.79 

25.45 

21.40 


11.00 

23.56 

19.30 

27.40 

34.10 

23.70 

22.81 

21.20 

23.00 

20.95 

26.31 

22.63 

— 

11.50 

25.72 

23.90 


36.26 

21.81 

19.92 

18.21 

18.40 

23.30 

29.25 

25.00 


12.00 i 

25.85 

24.82 

31.00 j 


27.95 

25.00 

26.05 

27.18 

25.30 

31.58 

25.80 


12.50 

28.39 

26.00 ; 

31.80 

38.90 

29.13 

26.57 

27.80 

28.19 

26.92 

33.91 

27.50 


13.00 

33.18 

27.70 

33.86 

41.25 

41.85 

28.42 

29.80 

30.08 

28.34 

35.68 

29.40 


13.50 

31.61 

29.66 

36.44 

44.08 

33.95 

30.70 

30.25 

32.74 

30.90 

39.77 

32.10 

— 

14.00 













14.50 













15.00 

33.06 

31.86 

38.00 


36.80 

34.02 

35.20 

35.16 

34.79 

41.95 

35.27 


15.50 













16.00 

28.00 

28.30 

33.85 


33.38 

31.10 

30.12 

32.00 

31.72 

38.60 

32.56 
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Figure 7.- Ratio of pressure to blunt flat-plate pressure without separation. 

occurring only in transitional separated flows. The present tests were conducted at local 
plate Reynolds numbers of approximately 100 and, in most cases, exhibited a region of 
fairly constant pressure. It is doubtful that transition could be influencing these results 
at a Reynolds number this low. 

Interaction length. - Figure 7 shows the values chosen to represent the constant- 
pressure plateau region. For those tests not having a definite plateau, the pressure at 
the inflection point near the hinge line was used. ■ The interaction length is defined, for 
this report, as the chordwise distance over which the pressure rise to plateau takes place. 
This distance was determined from plots such as figure 7 and compared in figure 8 with 
a theoretical expression derived by Erdos and Pallone. (See ref. 35.) The expression 
presented in reference 35 gives the interaction length in terms of quantities at the begin- 
ning of interaction and has the form 



(4) 
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1 2 3 4 5 6 7 8 



Figure 8.- Agreement between measured and predicted interaction length. 

where f^l) is empirically determined in reference 35 to be 3.43 for laminar flow, and 
| ^ is a correction factor to account for nonlinear effects at higher Mach numbers. For 
initial conditions (conditions at the interaction point) characterized by a flat-plate solu- 
tion, equation (4) can be represented by 


X f 4 (1)C P>P 


where for laminar flow N = 2 (from ref. 35). The form factor H 0 was evaluated in 
reference 35 by using the expression given by Monaghan (ref. 36). For the comparison 
in figure 8, the direct effect of the pressure gradient near the nose on 6 0 was 
neglected and equation (5) was used. The correction factor was assumed to be 
constant and determined by the best fit to the bulk of the experimental data. There is 
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considerable scatter about the line of perfect agreement. Perhaps this agreement 
could be improved by properly accounting for the effect of the history of the flow on 
the boundary -layer properties. 

Comparison of the plateau pressure with theory. - Figure 9 shows the measured 
plateau pressure coefficients along with some of the other existing laminar separation 
data. The local Mach number and Reynolds number at the interaction point were calcu- 
lated for the present tests by assuming an equilibrium isentropic expansion from the 
model nose. Calculations were also made for a flow frozen at the nose, but these calcu- 
lations showed no appreciable change in Rx,o and M 0 from the equilibrium case. 

The equation for the theoretical curves in figure 9 has been developed previously 
from linearized theory by several authors and has been experimentally verified at low 
Mach numbers. The equation has the form (ref. 8) 




Figure 9.- Effect of Reynolds number on plateau pressure coefficient. 
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At low Mach numbers with moderate boundary -layer cooling and small pressure gradi- 
ents, equation (6), for a flat plate, can be approximated by 

c p,p = B(R x ,o) _1/4 (/3o) _1/2 (V 

The constants A and B have previously been evaluated empirically and the value 
for B suggested by Sterrett and Emery (ref. 14) to fit the data of reference 8 has been 
used in figure 9. Considerable disagreement with equation (7) is seen at the higher Mach 
numbers. The present data are at a low local Mach number and are in fair agreement 
with equation (7). An attempt was made to use the more exact form of equation (6) by 
computing a skin-friction coefficient by using the method of reference 37. The agree- 
ment with equation (6) offered no improvement over equation (7). 

The present series of tests were conducted at a significantly higher stagnation 
enthalpy than the other data in figure 9. The introduction of this new variable may 
account for some of the scatter in the present data. 


Calculated Results for Extent of Separation 


The method described here allows the extent of laminar separation to be calculated 
with simplified physical concepts. This method is for wedge-induced separation where 


the separation and reattachment 
points are free to adjust them- 
selves to changes in flow condi- 
tions and body geometry. Fig- 
ure 10 illustrates the idealized 
flow model assumed for this 
method. 

Assumptions. - Several 
assumptions are necessary in 
the development of this technique 
and these assumptions are: 

(1) The pressure rises 
from interaction to plateau and 
from plateau to reattachment 
can be represented by oblique 
shocks 

(2) The entropy gradient 
between the flow in the separated 



boundary layer and the inviscid 
flow can be neglected 


Figure 10.- Flow schematic and pressure distribution. 


17 


(3) The dividing streamline is straight 

(4) The flow is parallel to the flap after reattachment 

(5) The total pressure along the dividing streamline at reattachment is equal to 
the static pressure downstream of reattachment as proposed by Chapman, et al. (ref. 8) 

(6) Separation occurs at the midpoint of the pressure rise to plateau pressure 

(7) The plateau pressure coefficient is adequately predicted by equation (6) 

(8) The interaction length is adequately predicted by equation (5). 

Most of these assumptions are straightforward and have been discussed by other authors. 
However, some discussion is perhaps necessary on assumption (5). The assumption that 
the total pressure along the dividing streamline at reattachment is equal to the static 
pressure downstream of reattachment is now open to question (ref. 38); however, for the 
purpose of this paper, it is assumed that assumption (5) is valid. Chapman applied this 
idea to the case of zero initial shear layer thickness. More recently this work has been 
extended by Denison and Baum (ref. 9) and Kubota and Dewey (ref. 10) to include the non- 
similar growth of the free shear layer with finite initial thickness. It is these more 
recent results that are utilized in the present method to determine the reattachment 
pressure. Assumptions (6), (7), and (8) were examined experimentally and were dis- 
cussed under the experimental tests. 

For the purpose of discussion, the method can be divided into two parts: the region 
around the separation point and the region around the reattachment point. A detailed 
step-by-step procedure is presented in the appendix and frequent references to steps 
in this procedure are made. 

Separation region. - When the flow properties at an assumed interaction point (see 
fig. 10) are known, the plateau pressure coefficient can be calculated by using equation (6). 
When Cp^p and M 0 are known, a flow deflection angle can be found by using oblique 
shock relations. For the case of a real gas, an iteration procedure ^given by steps (0) 
to @) in the appendix^ is required to determine the flow -deflection angle. The interaction 
length Zi can be determined by equation (5). The separation point is then fixed by con- 
sidering the pressure rise from p^. to p c to be a mirror image of the pressure rise 
from p g to Pp. Under these conditions, the separation point is given by 

x s = x o + J ( 8 ) 

and the dividing stream leaves the plate surface at x s with an angle 6 S . 

Reattachment region. - The equation for P r /Pp ^shown in step ^5) of the appendix j 
is taken from reference 8 and shows the reattachment or final pressure as a function 
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of y, local Mach number, and the ratio of dividing streamline velocity to local free- 
stream velocity just prior to reattachment. The value of y and Mach number are known 
since the pressure between separation and reattachment is assumed to be a constant. 

The only remaining unknown in the equation for p r is u*. References 9 and 10 provide 
graphs showing the variation of u* with a distance parameter S* for three different 
initial (separation point) profiles. The parameter S* is evaluated, as in steps (0) 
and 0), for a particular flap angle 9 and assumed interaction point. Then u* can be 
taken directly from the graphs for any of the three initial profiles and p r can be eval- 
uated as in step ( 0 ). For all cases presented in this report, results from all three pro- 
files are shown. 

Once the reattachment pressure has been determined, the flow -deflection angle 6 r 
can be found from oblique shock relations. It is obvious from geometry that 

6g + = 0 (9) 

Therefore, any solution must satisfy this equation, and a trial-and-error procedure is 
necessary to determine the correct flap angle for the assumed interaction point. The 
procedure in the appendix illustrates this trial-and-error technique. 

Comparison and experiment .- Figure 11 shows a comparison of the calculated and 
experimental results for the extent of laminar separation. The results are better than 
might be expected from such a simplified approach to the problem. Since 0° flap tests 
were not available for all tests, the undisturbed, blunt-plate pressures used to determine 
flow properties at the edge of the boundary layer for the theoretical curves in figure 11 
were obtained from the measured pressure distributions extrapolated to the plate -flap 
junction. The pressure gradient effect on skin friction was ignored and Cf 0 was com- 
puted by the Monaghan T' method. The plateau pressure coefficient was computed from 
equation (6). The exponential profile results seem to agree best with experiment. 

All the calculated curves in figure 11 are for equilibrium flow with the exception 
of the two dashed curves shown in figures 11(a) and 11(c). These two cases were calcu- 
lated for a flow frozen at the nose of the body. The frozen cases are seen to predict 
larger separated regions than the equilibrium cases at the same flap angle. 

Since the results obtained from the calculation technique were in fair agreement 
with the experimental tests, several calculations at various enthalpy levels and free- 
stream unit Reynolds numbers were made. Two additional assumptions were necessary 
to compute these cases. First, the undisturbed, blunt-plate pressure was computed by 
the method of reference 30 by using equations (1) and (2). Second, the wall temperature 
was assumed to be constant at 300° K. The skin-friction coefficient and plateau pressure 
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Figure 12.- Effect of enthalpy on separation length. = 12.0; T w = 300° K. 


coefficient were computed as before. Figure 12 shows the results of these calculations. 
This figure shows a large effect of free -stream unit Reynolds number and a smaller 
effect of enthalpy on separation distance. 

Not shown are several cases computed at different free-stream Mach numbers 
(10 < Moo < 20). Increasing the Mach number decreased the chordwise extent of laminar 
separation. 

Correlation parameter for exten t of s eparation . - The calculation technique just 
described has no closed-form expression for separation length. It would be helpful to 
see the manner in which the separation length depends on such variables as Mach num- 
ber, Reynolds number, skin friction, and model geometry. An expression of this type 
may be derived from a relation taken from reference 8: 


Z ds _ 1 

So* |/ c f,o )fP o 

when 


5 o 


* 


ClXp 


( 10 ) 


( 11 ) 
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( 12 ) 


Equation (10) then reduces to 


c f,o “ 


V^x~o 



(13) 


where 



Referring to figure 10 and assuming that 1^/2 can be neglected compared with 
yields 


m = - x Q jcos 6 S 

_ (l - x 0 )sin 6 S 
tan (0 - 6 s j ~~ 


(14a) 

(14b) 


^ds = m + n = ( L - x o) 

Combining equations (13) and (15) yields 

L - x O o 

X ° k\[&> 


cos 5 S + 


sin 6 S 
tan(0 - 6 S ) 


(15) 


(16) 


where 


£ = cos 6 S 


sin 5 S 

+ tan^0 - 6 S ^ 


(17) 


The calculation technique (in the appendix) provides all the necessary information 
for evaluating the parameter on the right-hand side of equation (16) for a chosen inter- 
action point. Figure 13 shows a plot of equation (16) where 


0 = 



(18) 
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Figure 13.- Correlation parameter. Moo = 12.0; hj = 2000 Btu/lbm (4.66 Mj/kg). 



and (p is evaluated from the calculation technique (in the appendix) for several different 
free-stream Reynolds numbers based on plate length. 

The term <p is seen to be not only a function of (L - Xq^Xq but also a function 
of Roo L* This dependence is illustrated in figure 14 where tfi is plotted against R^ l 
at constant (l - x 0 )jx 0 . Curves are shown for several different enthalpy levels. For 
500 Btu/lbm S h t S 3500 Btu/lbm (1.16 MJ/kg S h t ^ 8.15 MJ/fcg), and 10 3 £ Rco,l « 10 6 , 
the calculated dependence of <fi on R^ l can approximated by 
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<t> * ^(RoOjL) 1 ^ 5 Xq Xq = Constant^ (19) 

The 5000 Btu/lbm (11.65 MJ/kg) curve appears to be slightly displaced. Combining 
equation (19) with equation (16) gives 

■#■(*«, L ) 1/5 » < 2 °> 

Figure 15 shows the variation of <p(^R 0 o > with (L - x o)/ x o- Equation (19) is 

an approximation, but for most practical purposes it may be assumed that the variation 
of 4 >(Roo,l) 1/5 wit h (L - x o)/ x o reduces to a single curve for 10** ^ R,*, l = 10® 
and 500 Btu/lbm S h(A 3500 Btu/lbm (1.16 MJAg = ht 2 8.15 MJ/kg). 




x 


0 


Figure 15.- Modified correlation parameter. M» = 12.0; h( = 2000 Btu/lbm (4.66 MJ/kg). 

Figure 16 shows the present experimental data plotted in terms of the correlation 
parameter ^Roo^l)*^* The calculated curve was obtained from the calculation tech- 
nique given in the appendix. Figure 16 also shows the sharp-plate data of Miller, et al. 
(ref. 39), Johnson (ref. 40), and Putnam (ref. 33). The dashed lines are faired through 
the experimental data. There appears to be an additional dependence on Mq. The 
trend with Reynolds number and Mach number is in agreement with linear separation 
theory and at constant Mq the experimental data correlate reasonably well even for 
various wall temperatures. 
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Figure 16.- Correlation of experimental data. 

CONCLUDING REMARKS 

A study has been made of real-gas wedge -induced laminar separation. A series 
of experimental tests were undertaken at a Mach number of 12 and Reynolds number per 
foot of 10^ to establish the validity of a previously determined, perfect-gas relation for 
the plateau pressure coefficient. A calculation technique has been devised to predict the 
extent of wedge -induced, laminar separation utilizing simplified concepts. 

The results of the study indicate that the previously determined relation for pla- 
teau pressure coefficient as a function of local Mach number and Reynolds number pro- 
vides reasonable agreement with the present experimental data if real-gas conditions 
are used to evaluate local Mach number M Q and Reynolds number Q . The proposed 
calculation technique for the extent of wedge-induced laminar separation agrees reason- 
ably well with the present, real-gas, experimental data. The results of calculations made 
by using this technique show that for laminar separation (a) increasing the Reynolds 
number per foot increases the chordwise extent of the separation region, (b) increasing 
enthalpy above 2000 Btu/lbm (4.66 MJ/kg) also increases the size of the separated region 
at the Reynolds numbers investigated, and (c) increasing free-stream Mach number 
decreases the extent of separation. A correlation parameter was developed which indi- 
cated a dependence of length of separation on local skin friction and Mach number and 
Reynolds number based on plate length. The parameter correlates with fair success 
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some previous, hypersonic separation data at several values of the ratio of wall tempera- 
ture to total temperature. The trend of the results from the calculation technique devel- 
oped here for the correlated length of separation agrees with linear separation theory 
and, for a blunted plate, agrees with the bulk of the present experimental data. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., June 26, 1967, 

129-01-03-15-23. 


27 



APPENDIX 


STEP-BY-STEP COMPUTATION PROCEDURE 

This procedure shows the calculations necessary to compute the extent of wedge- 
induced laminar separation by the method outlined in this report. The procedure is suit- 
able for real air in equilibrium. To start the calculations, it is assumed that columns (T) 
to (4) are given. The remaining columns are evaluated in turn. Figure 17 defines the 
pertinent quantities used in the calculations. 



Figure 17.- Flow schematic showing velocity components. 


Column 

Heading 

Definition 

® 

Poo 

free-stream pressure, lbf/ft 2 (N/m 2 ) 

© 

Moo 

free-stream Mach number 

© 

P t,2 

stagnation pressure, lbf/ft ^ (N/m^) on model nose 

© 

h t,2 

stagnation enthalpy, Btu/lbm (J/kg) on model nose 

© 

P 2 

plate surface pressure prior to separation, lbf/ft 2 (N/m 2 ) 
(calculated by method of ref. 30 for this paper) 
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Column 

Heading 

Definition 

© 

h 2 

local value of static enthalpy at edge of boundary layer, 
Btu/lbm ( J /kg) 

© 

t 2 

local value of static temperature at edge of boundary layer, °R 


p 2 

local value of density at edge of boundary layer, slugs/ft 2 
(kg/m 3 ) 

® 

a 2 

local value of speed of sound at edge of boundary layer, fps 
(m/s) ((6), (j)> (8), and (9) may be computed by isentropic 
expansion from p^. 2 to P 2 on thermodynamic charts, 
refs. 41 to 44^ 

© 

p 2 

local viscosity at edge of boundary layer from reference 42 
at pg and T 2 , lbf-sec/ft 2 (N-sec/m 3 ) 

© 

v 2 

local velocity at edge of boundary layer, ^2(ht ? 2 - 1 * 2 ) or 
|^2^(4) - (6)j (Use appropriate units.) 

@ 

x o 

distance from model nose to beginning of pressure rise to 
plateau, ft (m). A value of x Q is chosen and calculations 
yield correct flap angle for this x Q . 

© 

Rx, o 

Pn V o x o (8)0 0 

local Reynolds number, ■ or where 

p o @ 

zero subscript quantities are local values obtained in 
steps (5) to (fj) at the chosen x Q 

@ 

© 

Mq 

V n 0 

local Mach number at x n , — ^ or XT 

a O ® 

0O = l/ M o 2 - 1 = \j@ 2 - 1 

© 

c f,o 

local-skin-friction coefficient at x Q (calculated by 
Monaghan T' method for this paper) 

© 


Ip 0 V G 2 or 1(8) @ 2 , lbf/ft 2 (N/m2) 
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Column 

@ 

@ 

@ 

© 

@ 

® 

@ 


Heading 

C P,P 


P P 


Definition 

plateau pressure coefficient at x 0 , 



plateau pressure, Cp pq Q + p Q or (O) (T?) + ( 5 ), lbf/ft^ 
(N/m2) 


" Pr 


/2c 


1.58 


f,o 


0o 


or 


w p/ w o 


p o/ p p 


W r 


'I'e 


normal velocity ratio defined in figure 17. A value of this 
ratio is assumed. 

Pn " Pn / W D \ 

static enthalpy in plateau region, h n + 1 + — or 

1 ° j 

local value of compressibility factor in plateau region from 
thermodynamic charts at Pp and hp 

local value of temperature in plateau region from thermo- 
dynamic charts at p p and hp, °R 

local value of density in plateau region at pp and Tp, 


ZpRTp 


or 


-, slugs/ft 3 (kg/m3) 


density ratio across oblique shock, 


© 


in 


If 


, then 


Pp is correct. If (£5) * Q2(|, then a new value of Wp/w G 
^ must be assumed and steps (2j) to (2Ej) repeated until 


velocity component defined in figure 17, 


, fps (m/s) 


2 (2D - (6 


1 - 


2(hp - h Q ] 

- ( w p/ w °y 


or 


shock angle defined in figure 17, arc sin w 0 /V 0 or 
arc sin 


separation angle defined in figure 17, 

\p a - arc tan^tan >/s)( w p/ w ojj or (27) - arc tan 


j(tan@)@ 
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Column 

® 

© 

© 

© 

© 

© 

© 


© 

© 

© 

© 


Heading 


Definition 


Mr 


Mr 


L w,o 


l aw,o 


H r 


A s 

L 

L - x c 


local velocity at edge of boundary layer in plateau region as 

defined in figure 17, . , Wp , or , 

fps (m/s) Sm ( fc ' 5s ) 

ratio of specific heats evaluated at Pp and Tp 

local viscosity in plateau region from reference 42 at Pp 
and Tp, (lbf-sec/ft 2 ) (N-sec/m 2 ) 


local speed of sound in plateau region, ^ypZpTpR or 
/(©©>@)R, fps (m/s) 


V D 

local Mach number in plateau region, — or 

a P 

plate surface temperature at x 0 , °R 

adiabatic wall temperature at x 0 , °R 

form factor at x G , 5 o *^0 o . Reference 35 suggests that H 0 
be calculated by 2. 59 


1 + 1.20|^° -l) + 0.4(2^i£-l 

i O / V 


or 2.59 


1 + 1.20 - 1 + 0.4 ^ - 1 

\® / V© ) 




interaction length defined in figure 10, 3.43x 0 H 0 

(from ref. 35) or 3.43 ©© p^^ ~ ft (m) 

l , 

distance from model nose to separation point, x 0 + ^ 
© + |©, ft (rn) 

distance from model nose to hinge line, ft (m) 


or 


distance from hinge line to separation point, 


, ft (m) 
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Column Heading 


Definition 



6 



x ds 



wedge angle defined in figure 17. 9 must be assumed and the 

following calculations done to determine whether the assumed 
6 is correct for the chosen x 0 . 


length of dividing streamline, (L - 



distance parameter from reference 


x s ^cos 5 S 

(m) 


(L - x g ) sin 5 S 
tan(0 - 5 S ) 


9, 


O.UppVpMp \ 

r x s x ds 

yjo p 2 V 2 m 2 cb y 


or 


0 . 


x=x s 

l 

vx=0 



. The step size Ax is chosen for 


desired accuracy or the denominator may be evaluated 
graphically. 


© 

© 


P r/ P P 


velocity ratio along dividing streamline from reference 10 

pressure ratio across reattachment region (from ref. 8), 

y P 



I. p ~ 1 M 2 
2 M P 


u* 2 ) 


y p“ 1 




or 



P r 



reattachment pressure, ^5)(I|), lbf/ft 2 (N/m^) 
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Column 

Heading 

© 

w r /w p 

© 

*\ 

h r 

© 

z r 

© 

^ f 

0 

P r 

0) 

Pp/Pr^ 

0 

’/'r 

® 

6 r 

0 

0 c 


Definition 

normal velocity ratio defined in figure 17. A value of this 
ratio is assumed. 


The same procedure as used in steps 0) to (0) is repeated 
for steps 0) to (0) until @ is equal to (0) 


shock angle defined in figure 17, arc sin w D /Vp or 
arc sin 


reattachment angle defined in figure 17, 

(tan@)@ 

@ +@. 

If 0) = 0), this solution is correct. If (0) ± (0), then 
a new value of 9 in 0) must be assumed and steps 0) 
to (0) are repeated until 0) = 0). 




arc tan 


tan y/. 


or 


- arc tan 


computed value of wedge angle, 6 S + 6 r or 
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